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Zero- Temperature Casimir Fluctuations 
and the Limits of Force Microscope Sensitivity 

J. A. SidlesB 

University of Washington School of Medicine 
(October 16, 1997 ) 

It is predicted that in force microscopy the quantum fluctuations responsible for the Casimir 
force can be directly observed as temperature-independent force fluctuations having spectral density 
97r/(401n(4/e)) h \8k\, where h is Planck's constant and Sk is the observed change in spring constant 
as the microscope tip approaches a sample. For typical operating parameters the predicted force 
noise is of order 10 -18 Newton in one Hertz of bandwidth. The Second Law is respected via the 
fluctuation-dissipation theorem. For small tip-sample separations the cantilever damping is predicted 
to increase as temperature is reduced, a behavior that is reminiscent of the Kondo effect. 
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The Casimir force is the mean force between two ob- 
jects that is generated by quantum fluctuations ID. In 
a recent review article 0, Barton notes: "It is strange 
that for nearly half a century after Casimir no curios- 
ity has been displayed regarding the fluctuations [of the 
Casimir force] about the mean." This lack of curios- 
ity is understandable in view of the prevailing opinion, 
as summarized in that Casimir fluctuations are "far 
too small to detect with any traditionally contemplated 
Casimir-type apparatus." 

In this article we venture a contrary prediction — that 
Casimir force fluctuations are large enough to be directly 
detected by force microscopes, and that these fluctua- 
tions provide a fundamental limit to force microscope 
sensitivity which is stringent enough to be significant in 
practical applications |[l| . We obtain this prediction via 
a strategy advocated by Kupiszewska || : 

The standard macroscopic quantum theory 
for nonhomogenous media . . . refers to a med- 
ium described by a constant refractive index. 
Although useful for many applications, this 
approach, as well as all other approaches ne- 
glecting losses, is generally incorrect. It is 
well known that the dielectric function must 
satisfy the Kramers-Kronig relations, other- 
wise causality would be violated. According 
to Kramers-Kronig relations, the imaginary 
part of a realistic, frequency-dependent di- 
electric function must not vanish, and that 
implies the dissipation of radiation energy. 
Therefore, a complete theory will have to in- 
clude not only the field and the atoms, but 
also a system that absorbs energy, usually 
called a heat bath or reservoir. 

In implementing Kupiszewska's program, we will con- 
sider the two force microscope geometries shown in Fig. ^. 
Both geometries assume a spherical cantilever tip. We 
confine our attention to experiments conducted in vacuo 
at cryogenic temperatures |0] , because such experiments 



FIG. 1. Two force microscope geometries: at left, tip vi- 
bration normal to the sample plane; at right, tip vibration in 
the sample plane. 

have the force sensitivity required to directly observe the 
predicted Casimir fluctuations. 

Our discussion centers upon four parameters which 
change when the cantilever tip approaches the sample — 
our goal is to predict these changes. The four parameters 
are: (a) the resonant frequency u>q, (b) the spring con- 
stant k, (c) the cantilever resonant quality Q, and (d) the 
force noise spectral density Sf. 

We begin by describing how these parameters are mea- 
sured. In typical experiments pj] , the tip is brought 
near the sample and the Brownian motion x(t) of the 
cantilever tip is observed interferometrically. In well- 
designed experiments x(t) is dominated by thermal noise, 
such that the cantilever motion is in thermal equilibrium 
with k(x 2 } — kgT, where (x 2 ) is the mean square tip 
displacement, kg is Boltzman's constant, and T is the 
ambient temperature. In such experiments the autocor- 
relation of x(t) is exponential: 



(x{t)x(t + T)) t = (x 2 ) e -"<W( 2C ?) cos(w r), 



(1) 



where () t denotes a time average. The autocorrelation 
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thus determines the parameters (a; 2 ), loq, and Q. From 
them, it is routine practice to: (a) infer the spring con- 
stant via k = mujQ, with m the motional mass of the 
cantilever, (b) verify that thermodynamic equilibrium is 
respected by checking that k(x 2 } = k B T, and (c) calcu- 
late a Langevin force spectral density Sf via 
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The term mwo/Q in (g) is recognizably the damping co- 
efficient of a mechanical oscillator; (0) thus represents 
the fluctuation-dissipation theorem |Q as it applies to 
force microscope experiments. By this method, or incon- 
sequential variants thereof, the four parameters ujq, Q, k, 
and Sf are routinely measured in force microscopy. 

For tip vibration normal to the sample plane (as in 
Fig. [I] at left), the spring constant k decreases as the tip 
approaches the sample. Our main prediction, which is 
derived in the second half of this article, is that the ob- 
served change in spring constant 6k will be accompanied 
by an increase in force noise SS f according to 



5S f = 



9tt 



h (-Sk). 
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Then from (^|) and (|^) , the dynamical cantilever damping 
muJo/Q is predicted to increase according to 
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thus ensuring that the Casimir fluctuations respect ther- 
modynamic equilibrium, such that k(x 2 ) — k B T for 
all values of 5k. Note that these predictions involve 
only experimentally measured quantities and Planck's 
constant — there are no model-dependent parameters. 

Are the predicted force fluctuations large enough to 
observe directly? From (|J), we compute that a spring 
constant shift 5k = —2.6 x 10 _3 N/m will be associated 
with a temperature-independent Casimir force noise of 
10 _18 N/v / Hz. Spring constant shifts of this magnitude 
are commonly observed in force microscopy. Recent ex- 
periments have demonstrated force noise levels of order 
10~ 17 NA/Hz, and if force sensitivity continues to im- 
prove jlTj then it is reasonable to expect that the pre- 
dicted Casimir fluctuations will become a dominant noise 
mechanism in future force microscope experiments. 

To describe experiments in which tip vibration is in 
the plane of the sample — as in Fig. [l] at right — we in- 
troduce a length scale I defined such that an attractive 
Casimir force / exerted between tip and sample gener- 
ates a change in spring constant 5k = f/l. To calculate I 
explicitly, let <j>(z) be the modal eigenfunction of the can- 
tilever, with z the coordinate along the cantilever length 
L, normalized such that <p(L) = 1. Then I is given by 



r 1 = 



/ dx 


d(f>{x) 
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Typically l/L is of order unity. Letting h be the tip- 
to-sample separation distance, the predicted increase in 
force noise as the tip approaches the sample is 
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There is no minus sign in this equation, in contrast to (^), 
because the end-on geometry shown at right in Fig. |l| 
yields a positive 5k at close tip-sample separation. Then 
(||) and (|^) yield the predicted damping increase: 

^ = =^Ar*. CO 
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Assuming 5k is independent of temperature to leading 
order — which is a reasonable assumption for Casimir 
forces — the predicted cantilever damping varies inversely 
with temperature, according to (|J) and (Q). Such inverse 
relations are uncommon in physics but they are not un- 
known; the Kondo effect is an example. 

It remains only to derive (|^) and (|^) by the program of 
Kupiszewska. We will not hesitate to make brutal simpli- 
fying approximations along the way, with a view toward 
obtaining results in closed form. Furthermore, we will 
finesse various model-dependent parameters by showing 
that they appear in 5Sf and 5k in such a manner that 
the ratio SSf /(5k) is parameter-independent. By this 
strategy we can reasonably hope to obtain final results 
which have broader validity than the underlying model 
from which they derive. 

Following Kupiszewska ^ , we model individual atoms 
as independent harmonic oscillators. By an appropriate 
scaling of variables the interaction of a tip atom a with 
a sample atom b can be described by the Hamiltonian 



H =\^a(pl + ql) + \u b (pi 
- PabqaQb + (heat bath). 



ql) 



(8) 

Here u> a and ui b are atomic oscillator frequencies and (3 ab 
is a dipole coupling whose strength and spatial depen- 
dence are discussed later. The operators {q a ,p a , qb,Pb} 
obey the usual commutation relations [g^, Pj] = ih 6{j , and 
it will turn out that no other information need be spec- 
ified about their physical nature. We specify the heat 
bath as the unique independent oscillator (10) model of 
Ford et al. |g] that induces linear damping in the Heisen- 
berg picture equations of motion: 

q a + T a u a q a + uj 2 q a = F a (t) + (i ab q b (9) 
q b + T b uj b q b + ulq b = F b (t) + [3 ab q a - (10) 

Here F a (t) and F b (t) are operator- valued Langevin forces 
which originate in the heat bath, and {r o , T b } are damp- 
ing rates. We pause to note that our conventions for cor- 
relation and power spectral density are those of Balescu's 
textbook Pi: 
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C ab {t) = - (A(t)B(t + r) + B(t + r)A{t)) 



Sa{u) 



dr e WT C AA (T), 



with ( ) an expectation over heat bath variables. Then as 
shown by Ford et al. Q , the assumption of linear damp- 
ing in (p|-[lo|) uniquely determines the Langevin force au- 
tocorrelation to be 



Cf^(t) 



2ir 



dw ftu coth( 



ftui s 
-2k B T> 



(11) 



The resulting model of atomic fluctuations resembles Ku- 
piszewska's model of gauge field fluctuations quite closely. 
Kupiszewska's model integrates over matter fields to ob- 
tain equations in which only gauge fields appear. Our 
strategy is the opposite: we have integrated over the lon- 
gitudinal gauge fields which generate the atomic dipole 
coupling (g), such that (p|-|lf|) contain only matter fields. 
The two approaches are formally equivalent, because in 
the real world gauge fluctuations and matter fluctuations 
are inseparable, such that either can be regarded as the 
fundamental dynamical variable. 

We now wish to compute the pairwise Casimir force f a t 
between two atoms, the associated spring constant k a b, 
and the force spectral density Sf a6 (a;). The force is given 
in terms of the gradient of the dipole coupling by 

(fab)i = (Vi/?« 6 ) (q a q b ), (12) 
from which k a b and Sf ab (lu) follow immediately as 

(kob)ij = -V i(f ab)j (13) 

(s u H).. = (v^xv./u) (i4) 

where the spectral density S( qaqb \(u>) is given by 



S 



(?a<?i, 



(15) 



The only dynamical quantities that appear in (JI4- 
|l5| ) are the expectation {q a %) and the autocorrelation 
C(q a q b )(q a q b ) ( T ) ■ Our next step, therefore, is to com- 
pute these quantities to leading and next-to-leading order 
in {f3 a b, T a , Tf,}, which are assumed small compared to 
{uJayUb}- Physically, the atomic oscillators are assumed 
to be weakly coupled by the Casimir interaction and un- 
derdamped by the ambient heat bath. 

As described in j^] and , any desired correlation in- 
volving Qg (t) and <7b(t) can be explicitly computed by 
solving (p|-|l l[) in the Fourier domain. For (q a <lb) an un- 
complicated Fourier integration yields 



(q a qb) 



ft/3 a b 



h(3 ab U!gUJ b 

,2\2 



2(u) a + Lu b ) L2tt(cj 2 -^) 

x( r a(l-4+ ln (4))+(«- & ) 



(16) 



plus 0(/3 3 ,/3F 2 ) terms. The zero of the denominator in 
( |l6| ) for ui a = Lot cancels against a zero of the numera- 
tor, such that (q a qb) is finite for all nonzero u a and w b . 
Setting r a = T b = in (|l6|), we recover the same expres- 
sion for {q a qb) as is obtained from the ground state of the 
Hamiltonian with the heat bath turned off. 

These results establish that heat bath damping does 
not alter the interatomic Casimir force in leading order, 
as is physically reasonable. 

With regard to the autocorrelation C( ?a g,,)( ?a g b ) (t) in 
(|i"5|), the Gaussian property of the Langevin forces, as 
discussed by Ford et al. || , allows it to be written as the 
product of two simpler autocorrelations: 

CW)^) (r) = C qaQa (r) C qbqb (r) + 0(/3 2 ab ). (17) 



In turn, C qaqa (r) is calculated by solving (|9|-|l l|) in the 
Fourier domain in a manner that precisely parallels the 
calculations of Li et al. Ml. The result is: 



„ , . T^o ft f°°, r Q UJ a U) COS(WT) 



nJa \luI-lo 2 + iT a u;\ 

ft u a r _r„i r i/2 f- \ 
- — e L ^ T " Z cos(w a r) 
2 uj a L 

+ llmg(^ a + ^)\r\)]+0((3l b ) 
^e- r ^ 2 cos(Lo a T)+0(T a ,p 2 ab ). 



0(f3 2 ab ) (18) 



(19) 
(20) 



By substituting a — * b we obtain C qbqb (r). Here u) a = 
(uj 2 — r^/4) 1 ' 2 is assumed real and positive, and g(z) is 
the exponential integral M defined by 



g[z) 



cos(f) 



dt 

t + z 



(21) 



The term Img((u) a + iT a /2)\r\) in (|g) can be proved to 
be monotonic in |t|, with initial value — tan _1 (r a /(2o) a )) 
and asymptotic value — r a (D a /(r 2 a^). This term there- 
fore describes squeezing of the quantum zero-point mo- 
tion by the heat bath damping — a phenomenon which is 
physically to be expected. 

The engineering import of ( p0| ) is that even at zero tem- 
perature, where classical oscillators exhibit zero noise, 
the zero-point motion of a lightly damped (r a <C u> a ) 
quantum oscillator carries ft/2 of noise power within a 
bandwidth r a centered on a carrier frequency u> a . 

Since each atom in our model is coupled to an inde- 
pendent heat bath, and is dynamically independent of 
adjacent atoms, we can obtain the total Casimir force 
and force noise by summing the atomic interactions pair- 
wise. For definiteness, we assume a Debye distribution of 
atomic frequencies {uj a ,u) b }, such that p(u>) — ?>u} 2 /u> D , 
where u>d is a Debye frequency — any other broad-band 
frequency distribution would yield similar results. Then 
we frequency-average (q a qb) from ([if]) and S( qaqb )(u>) 
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from (|T^), dl7|), and (|2(]) to obtain to leading order in 
{/3ah,r a ,r;,}: 



(<la,qb) 



du a dui b p{ui a )p{u b ) (q a qb) 



91n(4/e) hp ab 
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20uj d 



(22) 



(23) 



In (|23|) we have assumed u -C Wd, as appropriate for 
audio frequency force microscope experiments. 

It remains only to integrate over the tip and sample 
volumes shown in Fig. [l]. We specify that the tip and 
sample atoms have number density p a and pb respec- 
tively. The coupling f3 a b is assumed to have a dipole 
dependence /3 a b = ^/I r a6| 3 , with r a b the atomic separa- 
tion vector and k the strength of the dipole interaction. 
Substituting (|22]-|23]) into (|l2]dT4|) and integrating over 
the volume of a spherical tip of radius r that is separated 
by a gap h from a half-space sample, we obtain the total 
Casimir force, spring constant, and force noise as 



„ hp a pbK^ 

n 



31n(4/e)7r 2 



k*°* = -(n< 
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(24) 



(25) 



(26) 



Here {T)ij = <5y is the identity matrix, n is a unit vector 
normal to the sample surface, and (n <g> n)jj = n^nj is 
an outer product. With the neglect of O(hfr) terms, 
as is reasonable for close tip-sample approach, our main 
results (||) and (ph follow immediately from (p4]~|26]). 

Consistent with our stated goal of model independence, 
we need not specify numerical values for the parameters 
{p a ,Pb,K,ujE), T a , Tb, h, r}, because they do not appear in 
the final results. However, the dimensionless coefficients 
in (^) and (^|) are weakly sensitive to the functional form 
of the assumed Debye distribution of atomic frequencies. 
Thus, for example, the coefficient 97r/(40 ln(4/e)) ~ 1.83 
appearing in (j^) is best regarded as a coefficient of order 
unity, whose precise value will depend on the material 
properties and shape of the tip and sample. These coef- 
ficients are well suited for experimental determination. 

Arguably the two least realistic assumptions of our 
model are the assumed dynamical independence of ad- 
jacent atoms, and the coupling of each atom to an in- 
dependent heat bath possessed of an infinite number of 



degrees of freedom. The path to a more realistic model 
is clear but arduous. The heat bath model should be im- 
proved to describe realistic phonon and conduction band 
degrees of freedom, while taking into account the finite 
size of the tip, and electronic degrees of freedom in ad- 
jacent atoms should be realistically coupled. Both gauge 
fields and matter fields should be explicitly included in 
the Hamiltonian as in Kupiszewska's pioneering article 
H . The resulting field equations should be solved for re- 
alistic tip-sample geometries. Force fluctuations should 
be computed by the field-theoretic methods pioneered 
by Barton |9|JlO|] . Ideally, the results should explicitly re- 
spect the fluctuation-dissipation theorem and should be 
expressed in a simple and physically transparent form. 

Meeting these challenges will not be easy. Yet if the 
predictions of this article are experimentally confirmed, 
such that Casimir effects set the practical limits to force 
microscope sensitivity, then achieving a realistic under- 
standing of these effects will become a matter of urgent 
practical consequence, in particular to the biomedical re- 
search community [0. And if the predictions of this 
article are not confirmed, the question will be: why not? 

In either case, it is certain that Casimir effects will 
continue to engage and delight the imaginations of the 
physicists and engineers in coming decades. 
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